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ABSTRACT 

We develop an ansatz for expressing the free energy of the two dimensional g-states Potts model for 



C-~ i q > 4 near its first order phase transition point. We notice that for the moderate values of q < 15, the 

^^ ' energy profile at the phase transition is not expressible as a sum of gaussians. We discuss how this affects 

the traditional finite size analysis of this phase transition. In particular, the dominant length scale governing 



the finite size corrections turns out to be much ( ~ 6 times) larger than the largest correlation length in the 



^H. problem. 
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1. INTRODUCTION 

There has been a lot of recent interest in the numerical study of first order transitions. Though these 
systems are undeniably of great phenomenological importance, a purely analytical study has often been 
beyond the reach of current methods. The two-dimensional q-statcs Potts model is, in this arena, a very 
useful model: for, on the one hand, a lot is known about this model analytically [1- 6]; on the other, a 
number of numerical simulations [8- 18] have been performed on this model. 

Unfortunately, the finite size scaling behaviour of the numerical data are often not satisfactorily explained 
by the exact analytical results at hand. Small discrepancies also appear when the numerically measured 
surface tensions are compared to their theoretical estimates. 

In an attempt to clarify these issues, we tried to condense the analytically known results into an ansatz 
which might correctly incorporate the leading singularities of the free-energy near the transition point as 
function of q and /3, the inverse temperature. This allowed us to investigate all moments of the energy 
distribution function at the transition point. As expected from our previous analysis of a high q expansion 
in this model [19,20] and contrary to popular expectations, we discovered that the higher moments are very 
large at moderate values oi q ^ 10. As a result, the energy distribution cannot be satisfactorily expressed 
as a sum of gaussian distributions; and the usual methods of finite size analysis may have to be revised. 

In addition to this, the finite volume corrections appear to follow a simple scaling relation, where, 
however, the relevant length scale turns out to be much larger than the correlation lengths in the model. In 
fact only the simulations performed at a q value as large as 20 seem to have been in the truly asymptotic 
regime [16]. We believe that these non-negligible finite size effects (and proper accounting of the large energy 
cumulants) may have caused the apparent discrepancies in the numerical analyses. 

The organization of this paper is as follows. In section 2, we summarize the analytically known properties 
of the model. In the next section, we show how numerical methods allow to isolate finite volume effects 
directly in the free energy. We develop our ansatz in section 4, and use it to predict a set of scaling relations. 
We apply this ansatz to study the available numerical data in section 5. We end with a short discussion of 
the main results of our paper in section 6. 

2. KNOWN PROPERTIES OF THE MODEL AND PRELIMINARIES 

The g-states 2-D Potts model is defined by the Hamiltonian 

(y) 
where the <Ti variables, attached to the V = L^ sites « of a square lattice, take q distinct values; the 
summation extends to all the 2V pairs {ij) of neighbouring sites. We shall refer to V as the lattice volume. 
Many properties of this model have been known for a long time and can be found in the review by Wu[3]. 
Throughout this paper, we shall call 

F^ lim I — In z] (2.2) 

V ►oo V V^ / 



the free energy of the system, where Z is its partition function. The model has a temperature driven phase 
transition which occurs at a value of the inverse temperature /3 given by 

/3t = ln[Vg+l] (2.3) 

This value separates the region where the system is in either one of the q possible ordered phases (/3 > /3i) 
from the disordered region (j5 < /3t). When necessary, the quantities in an ordered or in the disordered phase 
will be labeled by a subscript (p — o oi d, respectively. 

The transition is second order for q < 4, where the critical properties are described by known, q 
dependent, indices. In particular, the specific heat and correlation length exponents a and ly are the same 
in all phases and are given by 

a{q) = 2{l-2u)/i3il-u)) (2.4) 

i^iq)^i2-u)/mi~u)), (2.5) 

where cos (§«) ~ \/q/'^- Note in particular that at q ~ 4, a and i^ coincide: 

a(4) = j/(4) = 2/3. (2.6) 

For each phase (p, we observe that although the correlation length £_,^ and the specific heat C^p both diverge 
for q < 4 when (3 — f3t > as 

(2.7) 

their ratio remains finite at q = 4 in the neighbourhood of (3t- In the expression above, the dots stand for 
(unknown) less singular and regular terms at /? = /?(. 

The phase transition a,t (3 — Pt is first order for q > A. The internal energies of the pure phases, E^, 
are exactly known at (3 — (3t [1]. The latent heat C ~ E^ — Eg, finite for any g > 4, vanishes as q — > 4^_; so 
that this quantity is continuous at the (second order) point q = A. Various expressions and properties of the 
Eip^s are given in Appendix A for completeness and further reference. 

In this first order region, the specific heats C<^ and correlation lengths ^ip are finite at any value of /3. 
The difference Cd — Co is known at (3t'- 

Cd-Co^ 13^ C/^ , (2.8) 

which shows that it vanishes in the same way as C when q — > 4+. However, their sum is unknown. Exact 
expressions for various correlation lengths have been obtained in the recent years [5,6] and details are collected 
in Appendix A. Here we quote only those results which are relevant to us. According to these calculations, 
and to their interpretation [6,7] in terms of the pure phase ^^'s, the disordered phase correlation length ^d 
has the following divergent behaviour as q — > 4_|_: 

U^^^^, (2.9) 
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where x is given by 

X = exp [n'^/i2e)] 



6* = In 



(2.10) 



The behaviour (2.9) was anticipated by Black and Emery a long time ago using renormalization group 
techniques [2]. Only very recently numerical simulations were able to produce data in agreement with this 
result [18]. 

When q — > 4+, 26 tends to as (g — 4)^/^ so that x diverges quite strongly. On the other hand, as q 
increases away from 4, 6 increases very slowly, growing only logarithmically as q — > oo, so that x stays very 
large over a very broad range of g-values. As a consequence, the approximation (2.9) of ^d remains valid 
over an equally large domain [7] (see Appendix A for the size of the corrections to (2.9)). For this reason, 
the variable x happens to be relevant to our forthcoming analysis of the phase transition above q — 4:. We 
note in particular that it also controls the leading behaviour of the latent heat, which vanishes as x^^'"^ at 
(7 = 4. 

So, the correlation length dX (3 — (3t in the disordered phase becomes infinite as q — > 4+ as it does at 
q = A when /? — > /3j_. Less is known about ^o, but it has been conjectured by Borgs and Janke [7] that 

Co = U/2 (2.11) 

which is consistent with numerical estimates of ^o at q — 10 [9,11]. 

If S,ip diverges at q — > 4+, it is not unexpected that the specific heat does so as well, since the latter is 
related to the integral of the energy-energy correlation function over the whole lattice: 

Here ^o is a fixed link, E^ (0 or 1) is the energy carried by the link t, l)^ denotes the average inside the 
phase (/3; and the thermodynamical limit V — > oo is taken before the limit j3 — > [3t. As it generically does 
at a 2"*^ order point, the sum in (2.12) most probably diverges along with the correlation length, the major 
contribution coming from the region of large ^o ~ ^ separation. Since S,d is known to diverge as q — > 4+ , Cd 
should also diverge. Further, since Cd — Co is bound to vanish in the same limit, Co also must go to infinity, 
consistent with the Borgs- Janke conjecture (Eq. 2.11). 

We continue this review of known properties by writing down the consequence of duality for the free 
energy F(/3) of this model. At /? > /3t (respectively /3 < /3t), F(/3) coincides with the ordered (disordered) 
free energy Fo (Fd). Duality relates F{P) to F (Pj , where /3 is related to /3 by 

(e'' ~l){c^-l)=q, (2.13) 



and wc have 



F, 



(^) 



In 



cxp 



(^ 



V9 



Fo(/?)-ln 



exp(/3) - 1 



V9 



(2.14) 



It is convenient to parametrize the distance to j3 — Pt (the transition point, where /3 = /3) by defining a 
variable h through 






(2.15) 



As (i — > [3t, b vanishes as (l + 1/-^^) (/? — /3t), the interchange fi < — > (3 is just 6 < — > —6, and duahty is the 
statement that the free energies of the model can be written 



Fo{l3) ^ Ft + b + g{b) 



FdiP) ^ Ft + b + gi-b) 



(2.16) 



Here Ft is Fg (/3t) — Fd (Pt), b is defined by (2.15), and g{b), which vanishes at & = 0, contains all the essential 
information about F{f3). It follows the n energy cumulants at f3 — f3t 



pi"-) 



d" 



^ - d/3"^^^^^ 



(2.17) 



/9=A 



once expressed in terms of the derivatives g^P'{0), p < n, oi g{b) have the same form for ip = o and d, up 
to the change g^^^O) < — > (— )''.g'^P^(0). Thus, the combination Fq — i^)"F^^' is completely determined by 
the knowledge of the cumulants of order lower than n. Equation (2.8) and the relation 



E, 



Eo = -2{l + iVg) 



(2.18) 



are two well known examples of this general statement of duality. 

Up to now we only considered properties of the model in the thermodynamical limit; in particular the 
limit V — > oo was always taken before the limit (3 — > Pt, which allowed us to properly define F{[3). We 
now turn to the properties of the system in a finite box with periodic boundary conditions, and recall the 
result obtained in Rcf. [4] for the corresponding finite volume partition function Zv{P): 



Zv{(3) = q cxp [VFo{/3)] + exp [VF^m + R{L, /3) 



(2.19) 



The first two terms represent the q ordered and the disordered phase contributions to the total partition 
function respectively, and the essential statement is that the rest R{L, (3) is bounded by the inequality 



\R{L,(3)\ < ci X q-"^"^ X exp(yF(/3)) 



(2.20) 



where ci and C2 > are constants. Similar statements hold for the logarithmic derivatives of Zv{P), that 
is for the cumulants of the energy distribution. We refer to R{L^ [3) as a "true" finite size effect, and to the 
first two terms in Eq. (2.19) as the asymptotic ordered and disordered contributions. 



We also recall that 

F(/3) = Foil3) for /3 > A 

F{I3) = Fd{p) for /3 < A , 
so that F — Max (i^o, i^c(). Strictly speaking, the result holds for q "large enough". We will ignore this 
caveat, and consider that (2.20) is valid for q> A. 

Eqs. (2.19,2.20) are very interesting since they control the size of any deviation of the partition function 
at finite size from the familiar two component representation Z = qZg + Zd- For this reason, we are going 
to demonstrate in the next section how F{I3) — i^t, close io (3 — [3t can, in principle, be extracted from very 
precise numerical data on the energy distribution, supplemented by the duality property (2.16). 



3. NON ASYMPTOTIC PURE PHASE FREE ENERGIES FROM NUMERICAL DATA 

By definition, the partition function is given by 

ZviP) = Y. ^v{E)eM~PVE] (3.1) 

E 

where VE is the energy of a configuration, and ily(ii^) (which is independent of /3), the number of con- 
figurations with that energy in the given volume. The current numerical simulations measure the energy 
probability distribution at some value of (3, say I3q: 

Pv,f}„{E) = nv{E)exp [-PoVE] jZy (M . (3.2) 

It follows from the two above equations that 

zvm 



Zv (/3o) 



J2Pvj,„{E)e-''^(^-M , (3.3) 



which yields the /3 dependence of ZyiP) up to a constant factor Zy (Po), and thus, an effective free energy 
via 

F^«'(i, /3) ^ F-'' (L, /3o) - -^ In [Zv(/3)/Zy (/3o)] . (3.4) 

This is true for any lattice system, but, in such generality, is not very useful for finite size studies because 
the V dependence of the constant (i.e. /3-independent) F'^^ {L,f3o) remains unknown. 

It becomes more interesting in the case of a system which has a first order transition at (3 = [3t, with, 
as in the Potts model, a duality property relating the q ordered phases to the disordered phase coexisting 
around [3t in a finite volume. Consider the "asymptotic" form Zas(/3) given by (2.19) when the remainder 
R{L,[3) is dropped and Fq^ Fd are represented in the form (2.16). Then we construct the quantity 

m - 4^.-M-W) ^ '^exp[y,(5)]+exp[y,(-6)] ^3^^^ 

■^as \Pt) 9+1 



with b given by (2.15) as a function of /3. As explained above, the function X{b) can be measured from 

numerical simulations in the vicinity of 6 = (that is of /3 — [3t). The function g{b) then follows from the 

knowledge of X at 6 and —6: 

-qX{b)^X{-b)- 



9{b) = ^\n 



(3.6) 
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and corresponding F^^ and Ff^ can be reconstructed from Eq. 2.16. By construction, g{h) should be 
independent of the lattice size in the limit where the remainder R{L, (3), bounded by (2.20), can be neglected. 
In turn, any V dependence is a direct measure of "true" finite size effects, and these effects, as opposed to 
the situation of the energy distribution (see our discussion in the next sections), are under theoretical control 

[4]. 

Let us now show that this extraction of finite size effects in the pure phase free energies is indeed 
feasible. For this purpose, we use the numerical data of Ref. [12], for the g = 10 Potts model with linear 
lattice sizes L = 16, 20, 24, 36, 44 and 50. Details on the runs can be found in the reference. We just recall 
that a multi-spin coding was used, allowing for the simultaneous generation of 16 independent Monte Carlo 
histories and that at each lattice size, a very long run was made near the maximum of the specific heat. We 
use the data from these long runs and all the errors quoted below are computed from the observed dispersion 
about their mean of the 16 corresponding results. Because it turns out that the finite size effects are very 
small, a plot of g{b) itself as given by Eqs.(3.6, 3.5 and 3.3) does not show any L dependence. For clarity, 
we choose to show, as a function of /3, and for the various lattice sizes, the quantity 

Wo[L, 13] = Ff{L, P)-Ft+ Eo (/3 - pt) , (3.7) 

which differs from the measured effective free energy Fo{L,f3) — Ft by the known asymptotic linear term. 
Then the finite size effects clearly show up, as exhibited in Figs, (la) and (lb). On the largest lattices 
(L — 50 and 44), because of important statistical errors, the quantity Wq is available in a limited range of (3 
only. Below [3 ~ Pt, this is due to the large cancellations occurring in the argument of the log in Eq.(3.6). 
At large positive values oi V {(3 — f3t) ^ the errors become large again since we have to extrapolate data far 
away from the (3 values where they were measured [12]. Fig. (la) shows Wo{L,/3) for L = 50, 44 and 36. 
From this figure, we conclude that, within the present errors, there is no sizable L dependence for L > 36, 
in the accessible (3 range, as was observed in Ref. [12] for (E) „ . 

On the contrary, finite size effects are clearly exposed (see Fig. lb) on lattices of linear size L < 36. In 
particular, the slope of Wo{L,(3) at (3 — (3t is not zero showing that the ordered phase internal energy is 
not Eg (but it does tend rapidly to E^ as L increases). This finite size effect should not be confused with 
that observed for the location of the peak in the energy distribution: as will be emphasized later on, a peak 
displacement must be present even in the asymptotic regime (i.e. even when R{L, (3) is negligible). A closer 
look at Fig. lb also reveals that the curvature (specific heat) at (3t (as well as the higher derivatives with 
respect to (3) do exhibit L dependence. 



The precision reached in the simulation used here does not allow for a quantitative study of the (3 
and L dependence of the "non asymptotic" part, R(L,(3), of the partition function, at least in the absence 
of some model for these effects. What is demonstrated is the feasibility of measuring the non-asymptotic 
contributions to the Potts model free energies. We now turn to a theoretical study of the asymptotic ordered 
phase free energy _F'o(/3), which the above numerical determination could be compared with. 

4. A THEORETICAL ANSATZ WHICH ACCURATELY DESCRIBES THE POTTS MODEL 
FREE ENERGY 

We start from the observation made in Section 2 that the correlation length of the model tends to 
infinity when the point q = A, (3 — /3t(4) is reached along cither of the following two paths in the q, /? plane: 

(i) q = 4, (3 > A (4) (Temperature driven 2"^^ order transition) 

(ii) (3 — (3t{q), q — > 4+ (along the 1^' order line). 

In case (i) we know furthermore that in each of the two phases ip = (o, d) the ratio of the specific heat 
to the correlation length tends to a constant (Eq.(2.7)). 

For case (ii), we argued that since ^.^ diverges, it looks reasonable to believe that C^ also diverges. 

We now conjecture that in fact dp/^ip is finite also along path (ii), and formulate our ansatz as a 
generalization of this kind of continuity of C^/^^ at q — A, l3 — /3t(4). (To be specific, wc discuss the ordered 
phase ((3 — > A+)j using duality afterwards to determine the free energy in the other phase.) 

In situation (i) the known specific heat exponent implies that the most singular part of the second 
derivative of the free energy with respect to (3 is 

Fi^Hp)^A{P-Ptr'^\ (4.1) 

so that the corresponding (p + 2) derivative is 

F(^+2)(/3)^^(-)^^^^±|P(/3-A)-'/^-^ . (4.2) 

On the other hand, we use our knowledge of the correlation length exponent: 

^o ^ Ao (/? -/?,)"'/' , (4.3) 

to eliminate {f3 — f3t) and obtain the relation 

F(^+^)(/3) ^ Ai-r^^^^±^ (eo/A„)-+^ , (4.4) 

between the order {p + 2) cumulant of the energy distribution and the correlation length. 

Our ansatz consists in boldly continuing Eq.(4.4), from situation (i) to situation (ii), where we know 
the behaviour of £,o (Eqs.(2.11,2.9)) as a function of q at /3 = (3t- Hence we have 

F(^+^) (A) ^ A(-)^^fc^(Bx)i+3^/^ , (4.5) 



where x is defined in Eq.(2.10). The constant B includes the proportionality constants occurring in ^o as a 
function oi (3 or oi x. Let us add a few comments. 

a) All these relations express the leading behaviour as either (3 — > /3t a.t q — 4, or q — > 4+ at (3 — f3f 
In particular, subdominant terms in the basic equation (4.5) arise from less singular or regular terms in both 

(2) 

Fo and ^o at the second order point q — 4:, /3 = /3i(4), neglected above. Including the (known) corrections 
to the low (q — 4) behaviour of ^o alone would thus have been inconsistent. 

b) We recall that x is very large over a wide range of values of q > 4. Accordingly, all the energy 
cumulants Fo (Pt), which behave like x(^+*^p/^\ become extremely large if p > 0. If we set p — —1, 
Eq.(4.5) yields a vanishing contribution of O (x^^''^) to the internal energy Eg. Though this must be 
neglected in comparison with other contributions to Fo (/?) which do not vanish at q = 4, /? — > /3t(4), we 
note that the known result 

Eo = -{l + 1/V?) - C/2, (4.6) 

shows the presence of this term as the latent heat C vanishes as a;^^ " (see Appendix) in the limit q — > 4+. 
e) At the leading order in x we consider here, the energy cumulants F^^ in the disordered phase are 
given by the expression (4.5) with the (— )^ sign dropped. This follows from duality (see Eq.(2.16)) which 
implies 

F^P+2)^(_)PFiP+2)+o(FiP+i)) , (4.7) 

where the corrections are down by x^'"^ compared to the leading term. 

d) In Eq. (4.5), we call the critical amplitude. A, and the proportionality factor, B, "constants". What 
we actually mean is that they are smooth functions of q when compared with the violent behaviour of x{q). 
This is indeed what appears in the latent heat C (Appendix), which is .t^^" times a function of q which 
varies smoothly from tt at g = 4 to An/ In q as q — > oo. 

e) For consistency with duality of our ansatz (4.5), we have the same A and B in the ordered and 
disordered phases cumulants. That Aq = Ad is no surprise since duality applied directly to (4.1) (which is 
vahd at (/3 - (3t) — > 0+) yields 

Ff (/3) ^ A{Pt^P)-^'\ 

which in turn is valid for [fit — /?) — > 0+. That it is so also for the factor B may seem more intriguing since 
^o is used in deriving (4.5) and is different from S,d- What it actually implies is that the ratio S,d/£,o at g > 4, 
which is two according to (2.11), must also be two at g = 4 (i.e. it must be continuous at q — 4+). 

With the expression (4.5) of the derivative Fo {j3) at Pt, n > 2, and the knowledge of the first order 
derivative —Eq, we obtain for Fo{0) — Fo {fit)'- 



'\A 

Fo{P) - Fo (A) = -Eo (/3 - Pt) + ^^ 



/ , \ 4/3 '^ 

(l + (/? - A) {Bxf/^) - - - (/3 - /?,) {Bxf/^ 



(4.8) 



For later convenience in exploiting this expression, we rearrange it and rescale its variables. After setting 

V ^ {(i ~ pt) {Bxf^ 

(4.9) 

Ed - Ed{BxY/^/{iA) , 
and noting that, to the order of accuracy retained, Fd follows from Fq by changing v into —v and Eq into 
—Ed,{*) we obtain the essential result of this section 

S^ (Fo(/3) ^Ft) = --^-v [e, + 1) + ^(1 + «)4/3 (4.IO0) 

S^ (FdiP) ~F,)^-^-v {sd - 1) + J(l - v^^' (4.10d) 

These are the expressions we propose as ansatze for the free energies of the model in the first order 
region. According to the way we established them, we expect their domain of validity to be the region where 
Bx is large and v is finite. Under these conditions, we immediately note that 

(i) S (which is proportional to the correlation length) plays the role of a length scale, to be compared 
with the linear scale L of the finite lattice on which we compute the partition function. 

(ii) After rescaling the internal energy parameters Eg, Ed according to (4.9) the "asymptotic" partition 
function 

Z - q exp {L^Fo{(3)) + cxp {L^Fd{p)) 

depends on L and q only through the reduced size L/ S (apart for the multiplicity q). This statement is true 
up to residual dependencies of A and B upon q in (4.9) (see comment (d) above). 

We repeat that the expressions (4.10o,4.10d) are just ansatze which have been built from the assumed 
continuity of ratios of energy cumulants at q = 4 and above, a continuity which we have no theoretical 
argument for. Hence, we have to justify our approach by explicit comparison of its consequences with results 
of numerical simulations. This will be done in the next section. However, we already want to point out that 
Eqs. (4.IO0 and 4.10d) suffer from the absence of singularity at f3 — f3t (that is u = 0). Such a singularity 
is expected [21] or proven to exist [22] in the similar situation of a field driven first order transition at low 
temperature (Ising model below Tc in the limit of a vanishing magnetic field). Here we take the heuristic 
point of view that this singularity, although probably essential, is numerically mild, so that it does not affect 
the practical consequences we shall now discuss. A probably related weakness of our ansatz is its singularity 
at u = — 1 for Fo ( and v = +1 for Ed). Indeed (1 + u)"'/^ is the only analytical function at w = whis has 
the required Fo (Pt)- But an essential singularity at w = 0, with a vanishing contribution to all derivatives, 
may alter the analyticity structure at w < 0. 



there is a misprint in Ref. 19 where Ed should be changed into —Sd in Eq.(12). 
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5. APPLICATIONS OF OUR ANSATZ FOR F{f3) 

In order to exploit the consequences of the ansatz (4.10), we need to assign numerical values to the 
"constants" A and B at the value of q we are interested in. As explained in Ref. [19], we are able to do so 
by comparing the values of the second and third derivatives of Fo{l3) &i [3 — [3t, given by 

i^i^) (^^) ^ ABx, 

2 (5.1) 

to the analytical results of a large q expansion. A full account of the procedure, which includes suitable Pade 
resummations of series in 1/ Jq now extended to order 10, is given in Ref. [23], together with estimates of A 
and B for the q values of interest. 

5.1. Comparison of Fo{/3) with numerical data at (j = 10 

As exhibited in Section 3, a simulation measuring the energy distribution in a region of /3 close to Pt, 
supplemented by duality, gives access to an estimate, Fg^{L,l3), of the ordered free energy. This estimate 
should coincide with Fo{P) in the asymptotic limit defined by the criterion that R(L,P) can be neglected in 
Eq.(2.19). In Fig. (2), we thus show the difference 

Ao = Ff{L, (3) - F^^"-^*") (/3) . (5.2) 

The quantity F^^{L,P) is the Fo extracted up to its constant term from the q = 10 simulation of Ref. [12] 
using Eqs. (2.16, 3.6, 3.5, and 3.3) while F™^^''' is the prediction of Eq. (4.10o), where we used the prede- 
termined values [19] 

A = .193 

B = .386 at q=lO (5.3) 

We did not plot the data for L — A4 and 50, which are always compatible with those for i = 36, with larger 
errors (see Fig. la). Let us now comment on Fig. 2. 

The lower limit of the /3 range shown corresponds to the value —1 of the scaled variable v of Eq. (4.9), a 
value below which our ansatz becomes meaningless. The upper limit corresponds to v about 2 and (3 values 
above which the data at L > 36 become unreliable. The /? range can also be expressed in terms of the 
variable {(3 — f3t) V , a natural argument for the partition function: it extends from about -4 to about 8 for 
L = 36. In this interval where we can compare the data and the ansatz, the following observations hold: 

(i) Although clearly visible, the finite size effects remain tiny: the maximum absolute value of L^Aq is 
less than 0.02 for all the data shown. 

(ii) For |f I < 1, the difference between the L = 36 data and the ansatz is less than 1 standard deviation. 

(iii) As already mentioned, a quantitative study of the (L, (3) dependence is difficult to achieve. The 
slope aX [3 — (3t is, however, measurable and shows a departure from the asymptotic value —Eq — 1.66425. 
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Wc find respectively at L = 16, 20, 24 and 36 

10^ X {Eo{L) - Eo) = -9.0(8), -5.1(1.0), -1.2(1.0), and 1.6(3.0) . (5.4) 

This clearly indicates a "fast" decrease of the finite size effect, of the same order of magnitude as the one 
found in [12] for the total average energy aX (3 — [3t, and compatible with a behaviour exp (— L/io) with 
Lo ^ 7. Note that this value is intermediate between S,d ~ 10.6 which should be relevant to the description 
of interface effects in R(L,f3) [7], and ^o ~ ^d/'2 ~ 5.3 expected to control the finite size effect in the pure 
ordered phase itself. 

We conclude that our ansatz does provide an accurate description of the asymptotic free energies in a 
substantial range of /? values around (3 — [3t. It can be used as a reference model for a quantitative study 
of finite size effects measured in a high precision numerical study for which the moderate size range L < 40 
would be especially suited. 

5.2. Energy distributions and finite size effects 

Given the free energies Fo{f3) and F^if)), the corresponding energy distribution can be computed by 
inverse Laplace transform of Eq. (3.1). As before, we call "asymptotic" the distribution obtained when 
R{L,f3) is neglected in (2.19), that is when 

ZviP) o^ Z^^p) = q exp [V Fo(/3)] + exp [V F4f3)] . (5.5) 

The probability distribution at f3 ~ f3o is then given by 



P^\{E) ^N d/3 exp [V (/3 - p^) E] Z^%P). (5.6) 

'J f3—ioo 

In this formula, A^ is a normalization factor fixed for eaeh V and f3o by 

J2PvME) = 1, (5.7) 

E 

(3 is suitably chosen according to analyticity properties of Z{/3). Note that strictly speaking, Fo and Ed 
should be periodic functions of complex f3 with period AiirV since the allowed values of E in Eq. (3.1) on 
a finite lattice are of the form E = —k/{2V), k ~ 0,1, . . . , 2V. Correspondingly, the integration contour in 
(5.6) should extend over a finite range, e.g., from /? — 2iTTV to /? + 2iTTV. We neglect such refinements: any 
reasonable Z^^ used in computing P{E) should not depend on details of F(/3) so far away from the real axis. 
This will be the case in what follows. 

Equation (5.6) of course implies that quadratic approximations of E^ and Ed as functions of (/3 — (3t) 
lead to the celebrated two gaussian formula [24] for the probability distribution. Higher terms in Fo and Ed 
generate size dependent deformations of the two peak structure observed in the vicinity of the first order 
transition, which should not be misinterpreted as 'non-asymptotic' contributions in the partition function. 
This caveat does not hold [4] if the energy cumulants a,t (3 — (3t only are analyzed. 
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Because, as we argued in previous sections, the Potts model free energies receive very large non quadratic 
contributions, F^"' (/3i) ^ x^^"^ for a broad range of q values, this model offers a valuable laboratory for 
studying their effect in measured probability distributions. With this remark in mind we will now compare 
to numerical data the predictions of (5.6) when Fq^ Fd are given by our ansatz (4.10o,4.10d). 

5. 2.1. Comparison with numerical data 

For dcfiniteness, we choose (3o — Pt in (5.6). The integral to be computed is thus: 

Pyp^ (E) ^ N dvx exp 

'^v — ioo 

x<gexp<l-) -w(eo + l) + 7(l + w)'*/^ ^ + exp 



^1 ev 





-v{+ed-l) + lil-vr/^ 



(5.8) 



Here we use the variables S, v, Eq and Ed defined in Eq. (4.9) and also scale the energy E according to 

e = EiBxy^y{3A) . (5.9) 

A'^ is a new normalization factor. This form (5.8) shows explicitly that, apart from the multiplicity q of the 
ordered phase contribution, P{E) depends on q and L through [L/ S) alone when expressed as a function of 
e variables. 

At large (L/S)'^ one may certainly evaluate the integral (5.8) by a saddle point method. The question 
we are trying to address here is how large should L be for this to be true. Taking q = 10 for dcfiniteness, 
using the numerical values from (5.3) and the expressions (4.9) and (2.10), we find 



S, 



(2=10 



60.3 



(5.10) 



about six times larger than the largest correlation length ^d^. Since we want to compare with distributions 
measured at L < 50, it is clear that the saddle point method is completely inadequate for this purpose. We 
shall use it later only for obtaining simple expressions on effectively large lattices. For the time being we 
have to perform a numerical integration of (5.8). 

Up to trivial changes, contributions in the disordered and ordered phase need similar analyses: only the 
ordered one is considered here. We have to choose v and a convenient path to ±ioo. We first remark that 
the stationary point of the integrand is 



vs{E)^-l + {eo-e + l)- 



(5.11) 



a value for which the integrand is well defined for 



Eo - e + 1 > 



(5.12) 



In such cases, we take v ~ vs{E) and the path v ~ v + i\, X E [— oo, +oo]. The term |(1 + v)'^'^ insures an 
exponential fall off of the integrand at large |A|. 
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For the other case, £ > Co + 1, we fix i; at —1 and again choose v — v purely imaginary. 

Our rcsuhs for L = 16, 20, 24, 36, 44 and 50 {with no arbitrary adjustable parameters) are shown on 
Figs. (3a, 3b) as continuous curves, together with the data of [12]. (Note the logarithmic scale.) The overall 
agreement, which continues over nearly three decades for P{E), is very good, it becomes nearly excellent at 
the largest L values. In particular the non-gaussian character of the peak is very well reproduced. 

We already know from our study of Section 3 that non asymptotic contributions are present when L < 36. 
We see their manifestation in Fig. (3a). Since the curves have unit area, the discrepancies between theory 
and numerical data are concentrated in the largest (the ordered) peak region, and are emphasized by plotting 
P{E) in that region on a non- logarithmic scale (Fig. 4a): the predicted peak looks too narrow. However, 
this does not mean that, at such L values, the effective ordered specific heat, or Fo {L,l3t), is larger than 
its asymptotic value: in fact Fig. 2 shows that the converse is true. At least for L = 16 and 20, the effective 
broadening of the peak arises from large asymmetries ( Fo {L,(3t) larger at low L values), whereas Fo 
is actually smaller (the curves in Fig. 2 are concave around j3t). Hence, when seen in energy distributions, 
the nature of finite size effects can easily be misidentificd. At the lower right corner of Fig. (4a), our curves 
level up; this unphysical feature is due to the tail of the ordered (left) peak. It disappears exponentially in 
L^ and is due to the singularity at u = — 1 in (4.10o) (see below the evolution of the ordered contribution 
far away from the peak). 

At larger L values (Fig. 3b), the agreement is considerably improved; in particular the absolute prediction 
for the ordered peak at L = 50 is perfect (Fig. 4b). However, a new effect tends to become sizable: a small 
discrepancy appears in the dip region between the two peaks a.t L ~ 50. This is the region where interface 
effects are expected to emerge [8], eventually yielding a fiat plateau \n E ai (3 = Pt- We shall soon discuss 
the effective L dependence of the dip depth, of direct physical interest for interface tension determinations. 

In order to understand better what is significant and what is not in finite size effects observed in the 
energy distributions, we now examine what the expected shape of Pv{E) is when L becomes actually large 
compared with the scale S* ( ^ 6 ^^ at g = 10) introduced in Section 4. 

5. 2. 2. Effective L dependence and large (L/S) limit of energy distributions 

Amongst the quantities which have been extracted, and tentatively interpreted, from numerically deter- 
mined energy distributions, are the peak locations E^2^{L) (asymptotically Eo,d) and widths (asymptotically 

related to Co,d), and the quantity 

'PviE"''''')' 



2aSf-^ln 



1 
L 



(5.13) 



Pv (E'^'P) 

(asymptotically twice the ordcrcd-disordered interface tension [8]). For definiteness, we discuss these quan- 
tities &i (3 — (it- 

Let us start with a study of the integral (5.8) via a saddle point method at large (, ~ L/S. Since by 
duality one peak is trivially deduced from the knowledge of the other, we again focus on the ordered peak, 
the first contribution in the integrand of (5.8). The saddle point vs{E) of (5.11) is useful when it lies above 
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the singular point of the ansatz, v = —I, that is in the energy region 

Rl : -2<E<Eo + 3A/{Bxy/^. (5.14) 

Ri extends above Eo by a small fraction (about 12% at q = 10) of the latent heat. It is the situation of 
interest for evaluating the asyniptotics of the peak structure. It is straightforward to find the saddle point 
result for the ordered phase energy distribution inside i?i. We find, up to an energy independent factor, that 

n4 



pordorocl^^ G Rl) = [1 + So - s] CXp <^ -t 



^'^'^-'^'--{l + e.-e) 



(l + 0(l/£2)) , (5.15) 



4 

e^L/s 

where we recall that e — Eq = (E — Eo) {Bx)^'"^ /3A. Note that P is asymptotically the exponential of a 
polynomial of degree 4. The prefactor corrects it by a term l/£'^ smaller. In order to analyze the approach 
to asymptotics of this peak, we expand the quantity 

Po{e) ^ -^ In P^^^^^^^iE) (5.16) 

around its maximum io, consistently up to 0(1) in l/l"^. The result is 

So - Eo - 1/ (3^') (5.17) 



Po{s) = est. + ^5 ( 1 + 72 ) ('^ ^ ^°)^ - (^ - ^o)^ 

^ ^ (5.18) 



^ '■ 1^ ' ~^4 , ^/ 1 . ~n5 



and 



4^1 + ^j(^-^~°) +0[^-^ie-eo) 

Corrections of order l/£^ are generic to the saddle point method. The above expansion is specific to our 
ansatz (note the absence of 1/1'^ term in the cubic term, and of any term surviving at L = oo at order > 5 
in e — io)- It predicts the approach to asymptotics of the ordered phase peak location, of the specific heat 
(quadratic term) and of the asymmetry (cubic term), and is relevant to a comparison with recent numerical 
simulations at g = 10 and 20 where such quantities have been considered [17]. Since £ is not large in most 
cases we also compare the latter results with those of the numerical integration of (5.8), from which we 
compute for each q and L the peak location and curvature at the peak. These comparisons are shown in 
Figures 5a, 5b, 5c and 5d for E°^^^ and Figures 6a, 6b and 6c for Co- 
Fig. 5a shows -B^ax for g ~ 10 as a function of 1/L (instead of the natural 1/i^) in order to emphasize 
once for all how misleading the finite size effects in energy distributions can be. Indeed, the data points 
of [17] as well as the continuous curve (our ansatz) exhibit an essentially linear behaviour in 1/L over the 
whole range L — 12-100. In the absence of the known exact value (square symbol) at 1/L = 0, a linear 
fit to the data points would be reasonable, but would predict an asymptotic value wrong by many standard 
deviations! This is not so for the q = 20 case (Fig. 5b), which reflects the fact that L/S is indeed large 
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for this case. Hence the apparent 1/ L behaviour clearly is a finite size artefact, and should not be taken 
as a significant property of the distributions [10]. Figs. 5c and 5d show the same quantities plotted against 
1/i^, together with the first order correction (straight line) predicted by Eq.(5.17). The adequacy of this 
representation is clearly exposed, as well as the strong deviations from asymptotics observed in the data and 
well (although not perfectly) explained by our ansatz. 

We choose to present Co (as extracted from the curvature at the peak maximum) as a function of 
S"^ / L'^^ in Figs. (6a) and (6b) for q = 10 and 20 respectively. Similar conclusions can be drawn. Note the 
ranges of S'^ jL? in the graphs are much larger than the domain of applicability of the saddle point method: 
S"^ / L"^ <C 1! As an illustration of the scaling property associated with Eq. (5.8), we finally show, on the same 
plot for g = 10 and 20, the ratio Co{L)/Co{<x>) against S"^ / L"^ (Fig. 6c). We see that the data points fall 
close to a same curve and to the theoretical prediction, even in regions where the asymptotic expectation is 
very far away. 

Let us now turn to the energy region in between the two peaks, which according to [8] can be (and 
has been) used for a determination of the interface tension between ordered and disordered phases. Let us 
recall the argument. Kt (3 — (3t and for sufficiently large linear size L, the mixture of a disordered phase in 
a fraction aV of the volume separated by two walls from an ordered phase occupying the fraction (1 — a)V 
has a relative weight which is 

i) independent of a because aFo + (1 — a)Fd = Ft at /3 = /3(. 

ii) proportional to exp [—2<7 odL] because each wall costs a free energy aL. 

Thus, a plateau in energy is expected, for any E of the form 

E = aEo + (1 - a)Ed, (5.19) 

whose height with respect to the rest of the distribution is exp [—2(TodL\. 

The predicted plateau has indeed been seen in recent simulations of the Potts model, either by extending 
the linear size of a square lattice to many times the correlation length (q = 20 in [17]), or by using rectangular 
lattices, which favour the appearance of separation walls transverse to the large dimension (q = 7 in [15, 11]). 
Such simulations have led to determinations of aod in very good agreement with the predicted (Jod = 1/ i'^Cd) 
of Ref. [7]. 

It is quite remarkable that previous numerical attempts on square lattices with sizes L < 100, not large 
enough at q < 10 for the plateau to develop, already led to satisfactory (although not as good) estimates 
of (Tod (Ref. [13] for q = 10 and Ref. [14] for q = 7). Indeed, as we saw above (Fig. 3a, 3b) our ansatz is 
able to reproduce the data for q = 10, L < 50 quite correctly, although it totally ignores contributions from 
mixed phase configurations. We illustrate this fact in Fig. 7a where we show our result and the data of [17] 
at g = 10 for the quantity 

1 ( / pmax pmax ~) 

2/o. = yln ^ ° " (5.20) 

^ I ^min J 
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as a function of 1/ L. The curve is computed from the mimerical integration leading to the energy distributions 
already shown. It levels out at low values of 1/ L (dotted part). This is because in the region 

^2: ^«+(^<^<^-(^' (5-21) 

the leading contribution to Po^d as L/ S becomes large, as obtained by a saddle point method, is proportional 



to 



exp 



L^\E-E,AI{Bxf"'] , (5.22) 



so that the right hand side of (5.20) grows as L. Only in this regime mixed phase contributions, of relative 
weight proportional to exp {—2(JodL) may show up, and this is supported by the data. What is puzzling 
is that for L < 40, our curve accounts for the trend of the data; moreover an extrapolation to 1/ L = of 
its nearly linear part would accidentally produce a reasonable result. The same remark applies to the case 
q = 7, illustrated in Fig. 7b (There the definition of fod is slightly different: we follow the definition of [14]). 
Our conclusion is that it may be misleading to identify fod, as defined by (5.13) or (5.20) or as in [14], to 
the interface tension cjod whenever a clear plateau is not seen in between the two peaks. 

6. CONCLUSIONS 

The main result of this paper can be summarized as that the ansatz for the free energy developed in 
section 4 is in qualitative and quantitative agreement with available numerical data. Our work indicates that, 
at least for the Potts model at moderate values of q, the energy distribution near a first order transition 
point cannot be accurately described as a sum of two gaussian contributions. This statement is distinct 
from the correct statement that, up to exponentially small (what we call non-asymptotic) corrections, the 
partition function can indeed be written as a sum of the partition function of the pure phases. 

These large energy cumulants of the asymptotic distributions call for more careful analysis of the 
numerical data. In particular, the finite size effects do not seem to be controlled by the correlation lengths 
of the model: rather, the relevant scale is significantly larger. This implies that numerical work should 
not assume that any system is sufficiently large: one can (e.g. using Eq. 3.6) and must specifically look 
for "true" finite size effects in the system (in the sense of Eq. (2.19) ). Especially, one must be careful 
about interpreting the depth of the minimum of the energy distribution as a measurement of the surface 
tension: true surface tension effects lead to a flat minimum which does appear on a sufficiently large lattice. 
In addition, deviations from predictions of the two-gaussian model (e.g. discrepancies in the position and 
curvature of the peaks in the energy distribution), should not, in itself be considered a finite volume effect: 
non-quadratic terms in the free energy expansion do lead to such shifts. 

The non-gaussian nature of the pure phase energy distributions near the transition point might seem 
unusual and pathological to the particular model at hand. In reality, our ansatz indicates that this is probably 
due to the influence of the nearby second order transition starting at q = 4. One may conjecture that strong 
deviations from the gaussian nature are not unusual in systems without a large surface tension between the 

16 



phases. Such might also be the case of phenonienologicaUy significant models like QCD, and the results of 
this paper might be relevant to any numerical simulation of these phase transitions. 

We thank J.M. Luck for his careful reading of our manuscript and for useful comments. 
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APPENDIX A 

For the q-statcs 2D Potts model on a square lattice, the internal energies Eg and Ed in the ordered and 
disordered phases respectively are given for g > 4 by [1] 

Ed + Eo = -2{l + l/^) (Al) 

C = Ed-Eo = 2{l + 1/Vg) tanh^P(0) {A.2) 

with 

2 cosh 6* = V? (-4.3) 

and 

oo 

P(6') = J]^ (tanh n6')2 . (AA) 

n=l 

Two expressions of P{9) [1] are of great interest for accurate numerical estimates of the latent heat C. 
One is very rapidly converging at large q: 

oo 

/p(^=l + 2^(-)™c-2™'V (A.5) 

?n— 1 

the other one is extremely accurate at small 9 

/ n \ 1/2 oo 
v/^PW=(y) ^ e-("+i) '^'Z^^^) (A6) 

\ ^ m=0 

In the paper, we chose to discuss all quantities in the low q^ A regime as functions of the variable x defined 
as 

a; = cxp [-7rV20] , (A.7) 

so that in this regime we have 

C^h{q)x-^''^[l + 2x-^ + 0{x-^)] (A.8) 

with 

Hq) - 2 (1 + 1/V^) tanh^ X {^^ (A.9) 

a smooth function compared to x. The corrections to the leading term in (A. 8) are extremely small, still less 
than 3% for q as large as 100. From these properties, we see that there is a very large domain of q values 
where large q as well as small {q — 4) expansions (A. 5) or (A. 6) yield very accurate estimates. 
Duality relates the specific heat difference to C: 

Ca-Co^Pl-^^o(x-^'^) . (AlO) 

Since we argue that Cd and Co are both ^ x at leading order, setting Cd = Co as implied by the final ansatz 
implies a relative error of order x^^''^ , not under control within our approach. 
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The interest of the x variable also appears in expressions for the correlation lengths [ 5, 6, 7]. 
According to the above references, the disordered phase correlation length S,d niay be accurately com- 
puted either at large q using 



1/e. = 4 Y. 

or at low (q — 4) using 



m— 1 



(-e-^) . , (me\ , , ^, ^, /cosh30/4\ 
-i '—smh tanh m0 + 2 log '— 



l/e. = 4^1n - 



1 + «;,; 



V2 cosh 



m+-)7rV0 



In particular the latter expression yields 



u 



8V2 






(All) 



(A.12) 



(A.13) 



subject to comments similar to those made on Eq. (A. 8). The accuracy of (A. 13) was already emphasized 
in Ref. [7]. 
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(Figure la) Graph of Wo[L, /3] (Eq. 3.7) at various large values of L, in an attempt to look for 
finite size effects a.t q = 10. Data from Ref. [12]. 

(Figure lb) Graph of Wo[L, (3] (Eq. 3.7) for small lattices a,t q — 10. Finite size effects clearly 
show up. Data from Ref. [12]. 

(Figure 2) Graph of Aq (Eq. 5.2), the difference of the measured free energy and the free 
energy predicted by our Ansatz. Data from Ref. [12]. 

(Figure 3a) Comparison of the energy distribution measured on small lattices against our 
prediction. Data from Ref. [12]. 

(Figure 3b) Comparison of the energy distribution measured on large lattices against our 
prediction. Data from Ref. [12]. 

(Figure 4a) Enlarged view of the peak region of 3a showing the difference between measured 
and predicted energy distributions. 

(Figure 4b) Enlarged view of the peak region of 3b showing the difference between measured 
and predicted energy distributions. 

(Figure 5a) Plot of the peak of the energy distribution versus 1/L at q = 10. A linear fit seems 
good, but predicts the wrong extrapolated value. Data from Ref. [17]. 

(Figure 5b) Plot of the peak of the energy distribution versus 1/L at q = 20. Data from 
Ref. [17]. 

(Figure 5c) Plot of the peak of the energy distribution versus 1/L^ at q = 10. Data from 
Ref. [17]. 

(Figure 5d) Plot of the peak of the energy distribution versus 1/L^ at q = 20. Data from 
Ref. [17]. 

(Figure 6a) Plot of the specific heat in the ordered phase at transition for q = 10 against the 
scaled lattice size, compared to our prediction. Here Co is extracted from the peak curvature. Data 
from Ref. [17]. 

(Figure 6b) Plot of the specific heat in the ordered phase at transition for q = 20 against the 
scaled lattice size, compared to our prediction. Here Co is extracted from the peak curvature. Data 
from Ref. [17]. 

(Figure 6c) Plot of the the scaled specific heat at transition compared to our prediction for 
both q = 10 and q = 20. Here Co is extracted from the peak curvature. Data from Ref. [17]. 

(Figure 7b) Comparison of 2fod defined in equation 5.20 with our prediction without any 
mixed phase contribution at q = 10. Data from Ref. [17]. 

(Figure 7a) Comparison of 2fod defined and measured in Ref. [14] with our prediction without 
any mixed phase contribution at q = 7. 
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